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Electron-scattering form factors for 6 Li in the ab initio symmetry-guided framework 
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We present an ab initio symmetry-adapted no-core shell-model description for 6 Li. We study 
the structure of the ground state of 6 Li and the impact of the symmetry-guided space selection 
on the charge density components for this state in momentum space, including the effect of higher 
shells. We accomplish this by investigating the electron scattering charge form factor for momentum 
transfers up to q ~ 4 fm _1 . We demonstrate that this symmetry-adapted framework can achieve 
significantly reduced dimensions for equivalent large shell-model spaces while retaining the accuracy 
of the form factor for any momentum transfer. These new results confirm the previous outcomes 
for selected spectroscopy observables in light nuclei, such as binding energies, excitation energies, 
electromagnetic moments, E2 and Ml reduced transition probabilities, as well as point-nucleon 
matter rms radii. 


INTRODUCTION 

The symmetry-adapted no-core shell model (SA- 
NCSM) p has been recently developed and designed to 
provide nuclear structure descriptions by using a new, 
symmetry-adapted and physically relevant many-particle 
basis. The model has been employed to unveil the emer¬ 
gence of a simple orderly pattern in nuclear dynamics, for 
the first time, in an ab initio framework (that is, from first 
principles), without a priori symmetry constraints. This 
highly structured formation is associated with an approx¬ 
imate symmetry in low-lying nuclear states that has been 
earlier suggested and linked to the symplectic Sp(3,M) 
group and its embedded SU(3) group [2lfT0]. The pat¬ 
tern favors low intrinsic spin together with large defor¬ 
mation and symplectic excitations thereof. This provides 
a strategy for determining the nature of bound states of 
nuclei in terms of a relatively small fraction of the pos¬ 
sible configurations. Consequently, we may extend the 
reach of ab initio approaches mm to explore ultra- 
large model spaces for a description of heavier nuclei and 
highly deformed structures together with the associated 
rotations. We have demonstrated that the SA-NCSM 
reduces the model space through a very structured se¬ 
lection, based on symmetry considerations, to physically 
relevant subspaces without compromising the accuracy 
of the ab initio NCSM approach pQ. 

In this paper, we focus on elastic (e, e') scattering 
charge form factors for the ground state of 6 Li and show 
that the SA-NCSM model with a symmetry-guided space 
selection provides a description of the form factors equiv¬ 
alent to the ones obtained in the corresponding complete 
space. This holds for any momentum transfer, from 
low q < 1 fm -1 through intermediate (up to 3 fm -1 ), 
and above (shown here up to q ~ 4 fm -1 ). While re¬ 


sults show that theoretical form factors are reasonably 
trending towards experiment, the 6 Li charge radius is 
not completely converged, so high-precision comparisons 
with experiment remain for future work. Nevertheless, 
the results presented here show, for the first time, that 
the calculated ground-state (gs) one-body charge density 
components in momentum space, including the contribu¬ 
tion from excitations to higher harmonic oscillator (HO) 
shells, is properly taken into account in selected spaces 
guided by Sp(3,M) and SU(3) symmetry considerations 
(similarly, for low-lying eigenstates of the gs rotational 
band). This, together with earlier SA-NCSM findings 
for observables such as binding energies, excitation en¬ 
ergies, electromagnetic moments, E 2 and Ml reduced 
transition probabilities, as well as point-nucleon matter 
rms radii for selected states [T], confirms the validity of 
the SA-NCSM concept. 

The significance of electron scattering form factors 
studies stems from their ability to provide a probe of the 
structure of the wavefunctions. For example, Ref. [2L 
studied inelastic scattering form factors and cross sec¬ 
tions to discern important spin flip components in 12 C 
wavefunctions that were sensitive to three-nucleon inter¬ 
actions. In this paper, we examine the longitudinal form 
factor (CO) for scattering off the ground state of 6 Li that 
is a Fourier transform of the ground-state charge den¬ 
sity. The CO form factors provide an indication on how 
well nuclear structure calculations reproduce the differ¬ 
ent lower- and higher-momentum transfer components of 
the nuclear charge density. This, in turn, can reveal im¬ 
portant underlying physics responsible for achieving con¬ 
vergence of the moments of the charge density starting 
with the rms radius. 

The charge form factors are calculated in the first-order 
plane-wave Born approximation. In all (e, e') calcula¬ 
tions presented here we use bare interactions, namely, 
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FIG. 1: iVmax = 12 model space (all circles, filled and unfilled) specified by the proton, neutron, and total intrinsic (S p , S n , S) 
spin values (vertical axis) across the Pauli-allowed deformation-related (A (i ) values (horizontal axis) for the 1 + ground state 
of 6 Li. The selected (6)12 SA-NCSM model space is shown by filled circles and it includes the complete space up through 
the 6 M2 subspace, while the 8M2, 10M2, and 12M2 subspaces are selected based on symmetry considerations favoring large 
deformation. The SA-NCSM results in a model space, for which all circles are filled up through the 12M2 subspace, coincide 
with the TV max = 12 NCSM results. 


Phil 2nil 


4hLl 


6 nn 


8M1 


io m 


12hLl 


- 



••••••••« 

oooooo»o«oi#• 

n i i i i i i i i i i i i i i i i i r 

oooooooooooto»oi#« 

niiiiiiiiiiiiiiiiiiiiiii— 

oooooooooooooooo«o«oi»»- 

- 




OOOOOOOOOOOOOo 

ooooooooooooooooo O 

OOOOOOOOOOOOOOOOOOOOOOOo- 

- 

• 



OOOOOOOOOOOOOo 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

OOOOOOOOOOOOOOOOOOOOOOOO- 

- 


... . 


o oo oo ooo o 

ooo ooo oo ooo o 

o oo oo ooooo o o ooo o - 

- 




0000000#0f0f## 

• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

- 




OOOOOOOOOOOOOI 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ooooooooooooooooooooooot- 

- 

• • • 



ooooooo»oto»»* 

• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

- 

• • • • 

. 


ooooooooooooo# 

Oooooooooooooooooo* 

ooooooooooooooooooooooor 





• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

• 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

- 




o oo oo ooo • 

ooo ooo oo ooo • 

o oo oo ooooo o o ooo • - 


-OOOJ — OO — coOcxj — O — O^ 

<o ct O' o' ct <o —~ <o c4 rt <o ct O cn <o ct in tj-\o oo O o~ ct -2 tj-VT <o ct in i> 'O oo o~ CD —~ <o O' ct -2 -st-VT o~ 'O rt in Tt tO 'so o\ oo o~rto'rt O ^ o'ci>n rf O cn \d o'> n ct ^foo' |V \oo\oo O o' rt ^ 


Deformations (A yu) 


the realistic nucleon-nucleon (NN) NNLO op t [22] and 
JISP16 [23] (with similar results obtained for N 3 LO [24]). 
The use of bare interactions, and not effective interac¬ 
tions in smaller model spaces, implies that operators used 
to calculate form factors does not have to be renormal¬ 
ized. In addition, charge form factors are calculated using 
the one-body charge density multipole operator, while 
contributions from two-body charge operators and/or 
relativistic corrections are not considered, as they are 
known to be negligible for charge form factors for mo¬ 
menta up to about q ~ 2 fm -1 [25] . Our calculated form 
factors have no center-of-mass (CM) contribution and are 
further adjusted to account for the finite proton size. 

SYMMETRY-GUIDED FRAMEWORK AND 

ELECTRON SCATTERING FORM FACTORS 

A detailed description of the ab initio symmetry- 
adapted no-core shell model (SA-NCSM) has been pre¬ 
sented, e.g., in Refs. (25| [27] . The SA-NCSM adopts 
the first-principle concept and is a no-core shell model 
(NCSM) carried forward in an SU(3)-coupled scheme [5j. 
The conventional NCSM m calculations are carried out 
in many-particle basis of Slater determinants (SD) built 
on HO single-particle states characterized by the M2 os¬ 
cillator frequency (or equivalently, the oscillator length 
b = ^/h/mQ). The model space is spanned by nuclear 
configurations of fixed parity, consistent with the Pauli 
principle, and truncated by a cutoff 7V max . The 7V max 
cutoff is defined as the maximum number of HO quanta 
allowed in a many-particle state above the minimum for 
a given nucleus. 

The many-particle basis states of the SA-NCSM for a 
given 7V max are constructed in the proton-neutron formal¬ 


ism and are labeled by the quantum numbers (A fi)nL of 
the SU(3 )(a m ) dSO(3)l group chain, together with pro¬ 
ton, neutron, and total intrinsic spins S p , S n , and S of 
the complementary SU(2) spin group. The label n dis¬ 
tinguishes multiple occurrences of the same L value in 
the parent irrep (A//). The orbital angular momentum 
L is coupled with S to the total angular momentum J 
with a projection Mj. Each basis state in this scheme 
is labeled schematically as |q N(\/jl)k,L; (S p S n )S] JMj ), 
where N is the total number of HO excitation quanta and 
7 denotes additional quantum numbers needed to dis¬ 
tinguish among configurations carrying the same TV (A /i) 
and ( S p S n )S labels. The organization of the model space 
allows the full space to be down-selected to the physically 
relevant subspace. 

The significance of the SU(3) group for a microscopic 
description of the nuclear dynamics can be seen from 
the fact that it is the symmetry group of the established 
Elliott model [3], and a subgroup of the Sp(3,M), the un¬ 
derpinning symmetry of the successful microscopic sym- 
plectic model 0. 

The charge form factors are calculated in the first- 
order plane-wave Born approximation. They are derived 
using the formalism and an extension of the computer 
code developed by Lee [28], described in detail in Ref. 
[25] . as well as using an SU(3)-based apparatus [29, 30 
for calculating charge and current density distributions 
in terms of the shell-model one-body density matrix el¬ 
ements (OBDMEs) and the single-particle matrix ele¬ 
ments of the associated electromagnetic operators. We 
calculate the OBDMEs using wavefunctions obtained in 
the ab initio SA-NCSM in complete TV max spaces or se- 
lected (A r m ax )Vn ax spaces. An {N^ ax )N max model space 
includes the complete basis up through N^ ax along with 
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FIG. 2: (Color online) Longitudinal CO electron scattering form factors F 2 (the Lo = 0 term in Eq. ([5]) with a procedure 
discussed in the text to produce a translationally invariant form factor) for the SA-NCSM 1 + ground state of 6 Li calculated 
in the complete iV max = 12 space (darker colors) and the SU(3)-selected spaces, (2)12, (4)12, (6)12, (8)12, and (10)12 (lighter 
colors), for Ml = 15 MeV or b — 1.66 fm (blue), Ml — 20 MeV or b = 1.44 fm (red), and HQ = 25 MeV or b = 1.29 fm (black) 
for (a) the bare JISP16 interaction, as well as for (b) the bare NNLO op t interaction. Experimental data are taken from Ref. 

m- 


selected (A fi) and (S p S n S) configurations beyond N^ ax 
up through 7V max (see Fig. 0 for a (6)12 model space). 

In the present analysis, we use the SA-NCSM with 
two realistic NN interactions, the bare JISP16 m and 
NNLOopt [22] potentials. The Coulomb interaction is 
added along with the NN interaction, together with a 
Lawson term for elimination of spurious center-of-mass 
excitations. We present results for 7V max = 12, as this 
model space is found sufficient to achieve convergence of 
the 6 Li gs energy - e.g., for INI = 20 MeV, it is within 
0.54 MeV of the extrapolated result of —31.49(6) MeV 
f3lff33] . Electron-scattering calculations are performed 
for a range of Ml = 15, 20, and 25 MeV and for several 
SU(3)-selected spaces, (2)12, (4)12, (6)12, (8)12, (10)12, 
together with the complete 7V max = 12 space. The result¬ 
ing wavefunctions, \aJMj) (where a distinguishes differ¬ 
ent eigenstates of given angular momentum J), are used 
to calculate lab-frame (or SD) OBDMEs, 

W J /ll{4iZiji;t z X ^ , n 2 l232\t z } ° || (1) 

where nlj label single-particle HO basis states and t z is 
either proton or neutron [a n ij-,t z is the annihilation SU(2) 
tensor operator, which destroys a proton or neutron in an 
nlj state]. These matrix elements are utilized to calcu¬ 
late longitudinal form factors for scattering from an ar¬ 
bitrary initial ( “i ”) eigenstate to an arbitrary final (“/”) 
eigenstate as a function of the three-momentum transfer 
q = |q|: 

f 2 l(q) = z^S~TT) ^ \< a f J f\\ M LoiQ)HJi)\ 2 , ( 2 ) 

Lo 

where the sum is restricted to \Ji — Jf\ < Lo < 
Ji + J f and M Lo (q) = Jj Lo (qr)Y^°(f)p(r)d 3 r = 


T,t=i( pvo tons)jL 0 (qri)Y^ 0 (fi) is the charge density mul- 
tipole operator given in position operators relative to the 
CM position operator. As the Ml 0 (q) is a one-body op¬ 
erator, its reduced matrix elements that enter in Eq. J2| 
can be expressed in terms of the OBDMEs of Eq. ([l]), 
provided the contribution of the CM component of the 
wavefunctions is properly removed. That is, the OB¬ 
DMEs are calculated for shell-model wavefunctions with 
lab coordinates and a CM component in the lowest 0s 
state, and hence, the CM-free F 2 (q) is calculated using 
the lab-frame OBDMEs of Eq. (HI) multiplied by an over- 

b 2 q2 * 

all factor of e 2 ^~, which removes the contribution of the 
CM component to the F 2 (q) [34] (see also [29]). In ad¬ 
dition, the calculated F 2 (q) form factors are adjusted to 
account for the finite proton size [35]. For the elastic 
electron scattering off the 6 Li ground state, the CO form 
factor, given by F 2 of Eq. Q with Lo =0, is calculated 
for Ji = 1, Jf = 1, and Jq = 0. 


RESULTS AND DISCUSSIONS 

Symmetry-guided form factors in the low- and 
intermediate-momentum transfer regime 

Longitudinal electron scattering form factors for the 
ground state of 6 Li are studied for the bare JISP16 and 
NNLOopt NN interactions up to V max = 12 spaces. An 
important result is that in all cases, (6)12 selected-space 
results are found to be almost identical to the V max = 12 
complete-space counterparts in low- and intermediate- 
momentum transfer regions (Fig. [2|, and even above 3 
fm -1 (not shown in the figure). This remains valid for 
various Ml values, as well as when different interactions 
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FIG. 3: (Color online) Longitudinal (70 electron scattering 
form factors E| (with a procedure discussed in the text to 
produce a translationally invariant form factor) for the SA- 
NCSM 1 + ground state of 6 Li calculated for HQ = 20 MeV or 
b = 1.44 fm and with the bare JISP16 interaction. The out¬ 
come for the SU(3)-selected spaces, (6)8 (red dots) and (6)12 
(blue dots), accurately reproduces the corresponding results 
for the complete iV max = 8 space (solid, red) and iV max = 12 
space (solid, blue), with larger-space iV max = 12 results lying 
slightly closer to experiment 36 . 


are employed (Figs. Hf and^). It also applies to ev¬ 
ery state of the 6 Li gs rotational band, as these states 
share the same SU(3) structure [l], but different total 
orbital momenta, and have very similar OBDMEs. This 
further confirms the validity of the symmetry-guided con¬ 
cept in the SA-NCSM. Indeed, while we have shown in 
Ref. [T| that the 7V max = 12 complete-space binding en¬ 
ergies, excitation energies, electromagnetic moments, E 2 
and Ml reduced transition probabilities, as well as point- 
nucleon matter rms radii are accurately reproduced in 
small selected spaces, the present results indicate that us¬ 
ing these selected spaces, that constitute only a fraction 
of the complete 7V max model space (about 1% for (6)12), 
reproduces, in addition, the 7V max = 12 complete-space 
form factor momentum dependence. In short, model- 
space selection, which is based on a straightforward pre¬ 
scription dictated by the Sp(3,M) and SU(3) symmetries, 
eliminates many-body basis states that are shown in this 
study to be also irrelevant for describing the charge dis¬ 
tribution for the 6 Li gs as revealed by the (70 form factor 
at low/intermediate momentum transfers and above. 

Deviations in the form factor as a result of the SU(3)- 
based selection of model spaces are found to decrease for 
higher Ml values (see Fig. [2j the higher the Ml value, 
the narrower the curve). This effect is more prominent 
for momenta q > 2 fm -1 . The outcome suggests that for 
high enough Ml values, results are almost independent 
from the model-space selection and, for Ml = 25 MeV, 
the (2)12 form factor already reproduces the 7V max = 12 


complete-space result. For low M1 values, larger 7\M ax 
spaces ((4)12 or (6)12) appear necessary pointing to a 
mixing of more deformation/spin configurations within 
these low-M2 spaces. However, while low values, h£t < 
15 MeV, are known to require larger model spaces to 
obtain convergence of the gs energy, such a mixing at 
the 4M2 and 6M2 subspaces is expected to decrease 
for 7V max > 12. In short, the SU(3)-based selection of 
the model space yields reasonably small deviations in the 
form factor, especially for q < 2 fm -1 and for Ml >15 
MeV. 

TABLE I: Binding energy (BE), excitation energies (E), elec¬ 
tric quadrupole (Q) and magnetic dipole ( fi ) moments, as well 
as point-nucleon proton (r p ) and matter (r m ) rms radii for 
the three lowest-lying T — 0 states in 6 Li, as calculated in 
the (6)12 SA-NCSM with the JISP16 NN interaction and for 
M2=20 MeV (taken from Ref. [l]) and compared to other 
ab initio approaches: the complete iV max = 12 model space 
[I] (or NCSM for JISP16 and M2=20 MeV), as well as Varia¬ 
tional Monte Carlo (VMC) and Green’s function Monte Carlo 
(GFMC) using the AV18 two-nucleon and Urbana IX three- 
nucleon interactions (energies taken from Ref. [41]; radii 
and electromagnetic moments taken from Ref. [42], without 
contributions from two-body currents). Experimental results 
(Expt.) taken from Ref. unless otherwise specified. 


SA-NCSM NCSM VMC GFMC Expt. 


1 + 

± gs 


BE [MeV] 

30.445 

30.951 

27.0(1) 

31.2(1) 31.99 

rms r p [fm] 

2.112 

2.125 

2.46(2) 


2.43 a 

rms r m [fm] 

2.106 

2.119 



2.35(3) 6 

Q [e fm 2 ] 

-0.08 

-0.064 

-0.33(18) 


-0.0818(17) 

g [^n] 

0.839 

0.838 

0.828(1) 

3+ 


0.822 

E [MeV] 

2.515 

2.526 

3.0(1) 

2.7(3) 

2.186 

rms r m [fm] 

2.044 

2.063 




Q [e fm 2 ] 

-3.88 

-3.965 




g \pn] 

1.866 

1.866 

2+ 



E [MeV] 

5.303 

5.066 

4.4(1) 

4.4(4) 

4.312 

rms rm [fm] 

2.18 

2.204 




Q [e fm 2 ] 

-2.279 

-2.318 




g [ gN] 

1.014 

0.97 





“Deduced from the 6 Li charge radius of 2.56(5) fm [36 
6 From Ref. 03] 

While results using NNLO op t lie slightly closer to ex¬ 
periment, both interactions show similar patterns with 
a small dependence on M2 (Fig. [2]). Furthermore, as 
one increases V max (e.g., from V max = 8 to V max = 12), 
SA-NCSM predictions are reasonably trending towards 
experiment, as illustrated for a (6)V max selected space 
and for the reasonable M2=20 MeV in Fig. [3j We note 
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FIG. 4: (Color online) Most dominant OBDMEs (with an absolute value > 0.02) labeled by ni,7i2(A/x) for a particle-hole 
(ni) 1 ^ 1 configuration, for the SA-NCSM 1 + ground state of 6 Li calculated in lab coordinates (with the Lawson term employed 
in the SA-NCSM calculations ensuring a 0s CM wavefunction component, ie., spurious CM excitations eliminated) and in the 
Amax — 12 complete space with the JISP16 bare interaction for hLl = 15 MeV (blue, left bars), 20 MeV (red, middle bars), and 
25 MeV (green, right bars). Error bars are defined by the range from the lowest value to the largest value of each OBDME 
over the set of SU(3)-selected spaces. 


that the A max = 12 results continue to deviate from the 
experimental data for intermediate momenta, especially 
for q > 2 fm -1 . Agreement with experiment may also 
depend on including contributions of three-body inter¬ 
actions in the SA-NCSM calculations and two-body op¬ 
erators in the F^. The significance of these contribu¬ 
tions has been shown in the framework of the Variational 
Monte Carlo (VMC) with the AV18 [32] two-nucleon and 
Urbana IX [40] three-nucleon interactions [41]. The low- 
HQ SA-NCSM F\ calculations using NNLO op t agree with 
the ones of the VMC using AV18/UIX (without contri¬ 
butions from two-body currents) for q < 2 fm -1 . The 
agreement might be a consequence of the fact that the 
NNLOopt is designed to minimize the contribution due 
to three-nucleon interactions (similarly, for JISP16). In 
order to gain additional insight into the similarities and 
differences among the ab initio results for 6 Li, we present 
in Table [I] the energies, electromagnetic moments, and 
point-nucleon rms radii for selected states in 6 Li, as calcu¬ 
lated in the present SA-NCSM approach with the JISP16 
and NNLOopt, and in other ab initio models, such as the 
VMC with AV18/UIX and the Green’s function Monte 
Carlo (GFMC) with AV18/UIX. The results presented 
in Table [I] show good correlations among the different 
models with, perhaps, the exception of the smaller rms 
radii obtained with JISP16 and the larger magnitude of 
the electric quadrupole moment obtained with the VMC. 
We note that the VMC with AV18/UIX has shown that 
two-body currents become significant for CO at momen¬ 


tum transfers of q > 2 fm 1 and are found necessary to 
achieve a close agreement with the experiment [4Tj . 

Important contributions to form factors 

One-body density for the ground state of 6 Li 

We study the most dominant OBDMEs for the ground 
state of 6 Li, as they are expected to provide impor¬ 
tant contributions to the form factor. By calculating 
the SU(3)-coupled OBDMEs (Fig. [4|, the largest ma¬ 
trix elements are found to belong to the njn^" 1 (A/i) = 
0 1 0 _1 (0 0) configuration (transitions within the s shell) 
followed by the l 1 l _1 (A/i) = (0 0) configuration (tran¬ 
sitions within the p shell), or An = \n\ — ri21 = Offl 
transitions. Typically, all the 0M2(0 0) contributions are 
important together with 2 h£t (A/i) = (2 0)/(0 2) and 
4 h£t (A/i) = (4 0)/(04), while there are smaller but non- 
negligible components for 6M2 (A /i) = (6 0)/(0 6) and 
0M2(2 2) (Fig. [ 4 ]) , followed by (A p) = (8 0)/(0 8) 
and 2M2(4 2) (not shown in the figure). 

The dominance of kh£t(k 0)/(0 fc), k = 2,4,..., in 
iV m ax = 12 complete-space OBDMEs (as shown in Fig. 
[4| can be recognized as another signature of the Sp(3, M) 
symmetry, as 2M2(2 0) single-particle excitations [and the 
conjugate (0 2)] are described by generators of Sp(3,M), 
while a stretched coupling of such excitations yields the 
multiples thereof [4M2(4 0), 6M2(60), etc.]. These, cou- 
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q(fm ) 


q(fm ) 


FIG. 5: (Color online) Most dominant SU(3) contributions to the translationally invariant Fl, denoted as F^ X ^ , for the 
longitudinal CO form factor, labeled by mri 2 (Xp) for a particle-hole (m) 1 ^ 1 configuration. The N max m 12 SA-NCSM 1 + 
ground state of 6 Li is calculated with the JISP16 bare interaction and for hQ = 15 MeV (left), 20 MeV (center), and 25 MeV 
(right). Note that the vertical axis scale is reduced by an order of magnitude from the top to the bottom panels. Results are 
very similar to the ones obtained with the NNLO op t bare interaction. 

*Components 6,2(4 0) and 7,3(4 0) lie almost on the top of the 4,4(2 2) and 5,5(2 2) components, respectively, for all q. 


pled to symplectic transitions of a particle two shells 
down, 2M1(0 2), can also yield 0M1(2 2), the result of 
(2 0) x (0 2), and 2Mf(4 2), the result of (4 0) x (0 2). 

In addition, we examine the dependence of the cal¬ 
culated OBDMEs on the symmetry-based space selec¬ 
tion (error bars in Fig. [4|. Specifically, the deviations 
are defined by the range from the lowest value to the 
largest value of each OBDME over the set of SU(3)- 
selected spaces, and are found to be reasonably small. 
Clearly, there is a very slight dependence on the model- 
space selection and on the value of HCt, with the exception 
of 2M4(2 0)/(0 2), which is about an order of magnitude 
smaller than the main (0 0) component. However, the 
comparatively larger uncertainties in 2Mf(2 0)/(0 2) are 
the reason, as also shown below, for the wider spread ob¬ 
served in Fig. [ 2 ] of the selected-space for momenta 
above q > 2 fm -1 . Moreover, OBDME amplitudes for 
2 1 0 -1 (2 0) and 4 1 2 _1 (2 0) (and conjugates) are found to 
decrease for smaller h£t values, eventually changing their 
sign. The observed opposite sign for small h£t has been 
suggested in Ref. m based on F 2 and charge densities 


of 6 Li and 12 C for hCl = 11 — 15 MeV. However, this 
effect appears to be independent of the type of the in¬ 
teractions employed, namely, the present study uses bare 
JISP16 and chiral interactions, while Lee-Suzuki effec¬ 
tive interactions for CD-Bonn and AV8' (plus a 3-body 
interaction) are explored in Ref. [25] . 


Momentum dependence of the form-factor SU(3) components 

Following Ref. HE EDI, where form factors are cal¬ 
culated in terms of SU(3)-coupled OBDMEs, we can 
study the SU(3) content of the corresponding electro¬ 
magnetic operators and the contribution of each (A p) 
term, F^ X ^, to the longitudinal form factor of the ground 
state of 6 Li as a function of the momentum transfer 
(Fig. [ 5 ]). The CM-free total longitudinal CO form fac¬ 
tor, F£, is given by the squared sum of all such terms, 

f 2 m = 

The results show that the largest contribution for all 
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yi(X /i) 




q (fm ) 


FIG. 6: (Color online) Same as Fig. [ 5 ] but for the NNLO op t 
bare interaction and for hQ = 20 MeV, with (solid) and with¬ 
out (dashed) removing the CM contribution (due to the Os 
CM component of the wavefunctions used to calculate the 
OBDMEs). See Fig. [ 5 ] for curve labeling. 


q values comes from the (A ji) = (0 0) (transitions within 
the s, p, sd, and pf shells), spreading to larger momenta 
for higher M2 (Fig. [5j top panels). As in the case of 
the OBDMEs discussed above, in addition to the strong 
0M2(0 0) contribution, the next important contribution 
comes from the 2M2(2 0) component (2 1 0 _1 and 3 1 1 —1 ), 
which peaks around 1-1.5 fm -1 (see Fig. [ 5 J middle pan¬ 
els, shown for a vertical axis scale an order of magnitude 



0.0 0.5 


1.0 1.5 2.0 2.5 3.0 

<7 (fm -1 ) 




FIG. 7: (Color online) Longitudinal CO translationally in¬ 
variant form factors, constructed from a subset of selected 
(A/i) OBDMEs, for the SA-NCSM 1 + gs of 6 Li calculated 
with the bare NNLO op t interaction (similarly for JISP16) for 
(a) HQ m 15 MeV, (b) hQ = 20 MeV, and (c) M2 — 25 
MeV, and compared to experiment (“+”) [36] . The (A/i) com¬ 
ponents [and (//A)] included are: all (grey solid), (0 0) (red 
dashed), (0 0) + (2 0) (red solid), (0 0) + (2 0) + (4 0) (green 
dot-dashed), (0 0) + (2 0) + (4 0) + (6 0) (blue dotted), and 
(00) + (20) + (40) + (60) + (80) (black long-dashed). Devia¬ 
tions due to the SU(3)-based space selection are indicated by 
the curve thickness. 


smaller than the one in the top panels). 

For intermediate-momentum transfers, for q > 2 fm -1 , 
Fl is predominantly influenced by 0M2(0 0), 2M2(2 0) 
(2 1 0 _1 and 3 1 ! -1 ), 4M2(40) (4 1 0 _1 and 5 1 ! -1 ), fol¬ 
lowed by 0M2(2 2) (2 1 2 -1 and 3 1 3 _1 ) (Fig. [5J bottom 






















panels). Compared to these contributions, the 6M4(6 0), 
2Ml(42), and 8Ml(8 0) components have a peak smaller 
in magnitude but located at slightly higher momenta, 
2.5-3 fm -1 , and become comparable in their contribu¬ 
tion around q 3 fm \ 

Furthermore, the changes associated with the SU(3)- 
based selection of the model space only appear to be 
significant for 2 1 0 -1 (2 0), and then only for low Ml (Fig. 
[5j widely spread curves). None of the other (A/i) con¬ 
tributions to the Fl are altered significantly by the 
SU(3)-based space reduction. In addition, for Ml =15 
MeV, a slight dependence on the space selection is ob¬ 
served for 3 1 1 _1 (2 0) as well as (but less importantly) for 
^(T^O), Pi- 1 (40), 2 1 2 _1 (22), and 6 1 0~ 1 (6 0), up to 
q < 2 fm -1 . However, for all Ml, the deviation observed 
for (4 0) is at least an order of magnitude smaller than 
that for (2 0). 

The effect of the CM Os component of the SA-NCSM 
wavefunctions on the form factor is illustrated in Fig. 
[6j As expected, the CM component suppresses the form 
factor due to its smearing of the translationally invariant 
charge density distribution. This effect has been demon¬ 
strated, for example, in 6 He [3J and 7 Li m- 

Finally, we consider form factors that are constructed 
of only several (A/i) contributions [together with their 
conjugates (/i A)]: starting with a form factor constructed 
of the (0 0) component only, and then consecutively 
adding the (2 0), (4 0), up to (8 0) components (Fig. [ 7 ]). 
Clearly, the (0 0) component makes up the predominant 
part of the form factor. It is interesting to note that, in 
this case, there is no dependence on the space selection 
for any Ml and for all q values. Also, for all Ml , the 
addition of the (2 0) component is found sufficient to re¬ 
produce the low-momentum regime of F\. Except for low 
Ml, the (0 0)+(2 0) form factor decreases at intermediate 
q values, while the consecutive addition of the (4 0), (6 0), 
and (8 0) components result first in an increase and then 
in a decrease of the intermediate-momentum F (Fig. 
[7] green dot-dashed, blue dotted, and black long-dashed 
curves, respectively). Those components are found to 
contribute the most to the form factor. In addition, for 
q > 2.5 fm -1 , including (2 2) to the F\ constructed of 
(0 0),(2 0), and (4 0), results in a slight increase of 
similarly, including (4 2) to the set of (0 0),(2 0), (4 0), 
(2 2), and (6 0) results in a slight decrease of F\. These 
(2 2) and (4 2) components slightly change the total 
and are found to be of a secondary importance. In short, 
for reasonable Ml values (> 15 MeV), the (2 0) compo¬ 
nent leads to a decrease in the intermediate-^ part of 
F\, bringing its value closer to the experimental data, 
while the (4 0), (2 2), and (8 0) are found to be foremost 
responsible to increase at intermediate momenta. 

We note that for smaller Ml values (< 15 MeV), re¬ 
sults (Fig. 00 are in agreement with the findings of 
Ref. [25] . Namely, the comparatively large (2 0)/(0 2) 
OBDME amplitudes in the wave functions, as discussed 


above, are found with the opposite sign to that needed 
to decrease F* and to reproduce the shape of the (e, 
e') form factors together with charge radii (the relation 

between the two observables can be seen from the low- 

2 

q expansion, F L {q 2 ) ~ 1 - (r 2 halge ) + ...). This has 
been clearly demonstrated in Fig. p] [the (2 0) panel], 
where the case of HQ = 15 MeV reveals a comparative 
large and positive (2 0)/(0 2) contribution for q > 1 fm -1 . 
The different behavior observed for low Ml is consistent 
with NCSM results for the 6 Li ground-state rms point- 
proton radius studied as a function of Ml and 7V max 
using the bare JISP16 NN interaction [31]. This study 
has revealed that for Ml < 15 MeV, the radius exhibits 
a larger dependence on Ml, while a steady increase with 
7V max (implying a decrease for F^) is observed only for 
hfl > 15 MeV. The importance of the (2 0)/(0 2) OB¬ 
DME amplitudes, their Ml-dependence and sign (known 
as The sign problem’ [25], not to be confused with the 
term, e.g., used in Monte Carlo approaches), merits ad¬ 
ditional investigation including their roles in other states 
and other nuclei [38] . 

CONCLUSIONS 

Longitudinal electron scattering form factors for the 
ground state of 6 Li were studied in the framework of 
the SA-NCSM for the bare JISP16 and NNLO opt NN 
interactions for a range of Ml = 15,20, and 25 MeV 
and for several SU(3)-selected spaces, (2)12, (4)12, (6)12, 
(8)12, (10)12, together with the complete 7V max = 12 
space. An important result is that in all cases, (6)12 
selected-space results are found to be almost identical to 
the 7V max = 12 complete-space counterparts for any mo¬ 
menta, shown here up to momentum transfer q ~ 4 fm -1 , 
while being reasonably close to experiment. This remains 
valid for various Ml values, as well as when different 
bare interactions are employed. Deviations in the form 
factor as a result of the SU(3)-based selection of model 
spaces are found to decrease for higher Ml. This effect is 
more prominent for momenta q > 2 fm -1 . However, for 
high enough Ml values, results are almost independent 
from the model-space selection and, for Ml = 25 MeV, 
the (2)12 form factor already reproduces the V max = 12 
complete-space result. 

The outcome shows that the largest contribution comes 
from the (A/i) = (0 0) OBDMEs and, for all q values, 
from the associated (0 0) contribution to the Fl, which 
makes the diagonal one-body density (within the s, p, 
sd , and pf shells) most important. In addition, the Fl 
for higher momenta, q > 1 fm -1 , is also influenced by 
2Ml(2 0) (2 1 0 _1 and 3 1 !- 1 ), 4M2(40) ^O” 1 andS 1 ^ 1 ), 
followed by 0Ml(2 2) (2 1 2 _1 and 3 1 3 _1 ) and also 6M2(6 0) 
(6 1 0 _1 and 7 1 ! -1 ). There is a very slight dependence on 
the model-space selection and as one varies the value of 
Ml, with the exception of the 2Ml(2 0)/(0 2) component. 
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However, the 2M)(2 0)/(0 2) OBDMEs are about an or¬ 
der of magnitude smaller than the those for the main 
(0 0) component. In addition, for all M2, only the (0 0) + 
(2 0)/(0 2) components are found sufficient to reproduce 
the low-momentum regime of F\. The (4 0), (2 2), and 
(8 0) components are the ones that are most responsible 
for larger F\ values at intermediate momenta. The pre¬ 
ponderance of 0M2(0 0), 2M2(2 0), ..., and 8M2(8 0) to¬ 
gether with 0M2(2 2) and 2M2(4 2) (and their conjugates) 
in the OBDMEs as well as the associated contribution to 
Fl can be recognized as another signature of the Sp(3, M) 
symmetry. 

In short, model-space selection based on Sp(3,M) and 
SU(3) symmetry considerations of the type we consider 
in the symmetry-guided concept of the SA-NCSM and 
that has been used to describe the low-lying structure 
of 1 2 3 4 5 6 7 8 Li in Ref. [I], properly treats, in addition, the 6 Li 
ground-state form factor for any momentum transfer 
(shown here up to g ^ 4 fm -1 ). The symmetry-adapted 
model spaces include the important excitations to higher 
HO shells as seen in their significant contributions 
at low- and intermediate-momentum transfers. The 
outcome further confirms the utility of the SA-NCSM 
concept for low-lying nuclear states. 
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